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ABSTRACT 

To determine the relative contributions of galactic and intracluster stars to the enrichment of the 
intracluster medium (ICM), we present X-ray surface brightness, temperature, and Fe abundance 
profiles for a set of twelve galaxy clusters* for which we have extensive optical photometry. Assuming 
a standard IMF and simple chemical evolution model scaled to match the present-day cluster early-type 
SN la rate, the stars in the brightest cluster galaxy (BCG) plus the intracluster stars (ICS) generate 
31^9^%, on average, of the observed ICM Fe within rsoo (~ 0.6 times r2oo, the virial radius). An 
alternate, two-component SN la model (including both prompt and delayed detonations) produces a 
similar BCG-HlCS contribution of 22^^%. Because the ICS typically contribute 80% of the BCG-hlCS 
Fe, we conclude that the ICS are significant, yet often neglected, contributors to the ICM Fe within 
r5oo- However, the BCG-flCS fall short of producing all the Fe, so metal loss from stars in other 
cluster galaxies must also contribute. By combining the enrichment from intracluster and galactic 
stars, we can account for all the observed Fe. These models require a galactic metal loss fraction 
(0.841q :[4) that, while large, is consistent with the metal mass not retained by galactic stars. The 
SN la rates, especially as a function of galaxy environment and redshift, remain a significant source 
of uncertainty in further constraining the metal loss fraction. For example, increasing the SN la rate 
by a factor of 1.8 — to just within the 2a uncertainty for present-day cluster early-type galaxies — 
allows the combined BCG -|- ICS -I- cluster galaxy model to generate all the ICM Fe with a much 
lower galactic metal loss fraction (^ 0.35). 

Subject headings: galaxies: clusters: general - X-rays: galaxies: clusters - galaxies: intergalactic 
medium - stars: supernovae: general 



1. INTRODUCTION 

The high measured Fe abundance of the intracluster 
medium (ICM) of galaxy clusters, ~ 0.3 "Zpe^Q (Edge 
& Stewart 1991), has been challenging to explain. The 
source of these metals is still unknown. Various enrich- 
ment mechanisms have been proposed, including models 
in which cluster galaxies lose enriched gas via galactic 
winds (De Young 1978), ram-pressure stripping (Gunn 
& Gott 1972), or gravitational interactions (Moore et 
al. 1996), but these have fallen short of producing the 
observed metals unless non-standard assumptions are 
made; for example, top-heavy initial mass functions 
(IMF) (Loewenstein & Mushotzky 1996; Gibson & Mat- 
teucci 1997), severe mass loss from galaxies (Renzini et 
al. 1993; Renzini 1997), or pre-enrichment from Popula- 
tion HI hypernovae (Loewenstein 2001). 

Such exotic explanations are intriguing, but we must 
first rule out more mundane explanations. One source of 
ICM metals that has only recently been considered is the 
diffuse, stellar component in clusters, i.e., the stars that 
lie outside the member galaxies (Domainko et al. 2004; 
Zaritsky et al. 2004; Lin & Mohr 2004). These stars must 
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be a source of metals for the ICM because they pollute 
in-situ; the key question is whether they are a significant, 
or perhaps even the dominant, source. 

Recent measurements of the intracluster stars (ICS) 
in a large sample of clusters demonstrate that this com- 
ponent contains roughly 30% of the stellar mass out to 
rsoo^ and an even larger percentage at the smaller radii 
typically probed by X-ray measurements (Gonzalez et al. 
2007). In previous work, Gonzalez et al. (2005) detect 
these stars in every cluster in their sample out to a radius 
of at least 300 kpc {h = 0.7). By stacking SDSS clus- 
ter images, Zibetti et al. (2005) detect the averaged ICS 
to 700 kpc, demonstrating the large extent of this com- 
ponent. Without any further calculation, one concludes 
that these stars could contribute ~ 30% of the ICM's Fe. 

Our first test of this scenario was promising, suggesting 
that intracluster stars could play a significant role in ICM 
enrichment (Zaritsky et al. 2004). This result was pre- 
liminary because we lacked spatially-resolved metallicity 
profiles to match apertures between the optical and X-ray 
measurements, and we had the minimum necessary data 
to make the optical/X-ray comparisons for just three 
clusters, of which only two had intracluster light mea- 
surements based on our high sensitivity methods. Be- 
cause intracluster stars are typically more concentrated 
than the cluster galaxies (Gonzalez et al. 2007), and clus- 
ters often have radial abundance gradients (De Grandi & 
Molendi 2001), our test must be re-done using matched 
apertures and homogeneous data. 

^ The radius within which the mass overdensity is 500 times the 
universal value and which corresponds to ~60% of the virial radius 
(Cole & Lacey 1996). 
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Wc have measured uniformly the properties of the in- 
trachistcr stchar components for a sample of 28 clusters 
(24 of which are from Gonzalez et al. (2005), while the 
remainder are from Kelson et al. (2008)). In this paper, 
we use XMM-Newton archival and guest observer data 
to measure the radial distribution of the X-ray emitting 
plasma and the Fc within that plasma for those 12 clus- 
ters with available X-ray data. We then determine the 
gas and Fe masses enclosed within the same radii used 
for the stellar mass measurements. 

The outline of this paper is as follows. In §2 we review 
our sample. In §3 wc discuss the data and our analysis. 
In §4 we present the X-ray spectral, morphological, and 
gas mass measurements. In §5 we discuss the compari- 
son of the observed Fe abundances with those predicted 
by a simple chemical evolution model assuming metal 
enrichment from the intracluster and galactic stars. Fi- 
nally, in §6 we summarize our conclusions. Through- 
out this paper, except where comparing to previous work 
and noted, we adopt the concordance cosmological model 
(Qa = 0.70, rirn = 0.30, and Ha = 70 km s'^ Mpc'^), 
which we call LCDM70 and dc;finc tx to be the radius 
of a sphere within which the mass density is X times 
the universal mean. All reported errors in this work are 
quoted at the la level. 

2. SAMPLE 

The clusters in this study are a subset of those for 
which we have intracluster light measurements (Gonza- 
lez et al. 2005; Kelson et al. 2008). First, we searched 
the XMM-Newton Science Archive (XSA) for European 
Photon Imaging Camera (EPIC) observations of any of 
these clusters. Observations of eight clusters arc avail- 
able, although one, Abell 3376, is not included here due 
to its highly disturbed X-ray morphology. Second, we 
obtained XMM-Newton Guest Observer (GO) observa- 
tions of five additional clusters. The selection criteria for 
the sample of twelve clusters presented here (see Table 
1) are a combination of those in our original intracluster 
stellar light work, those by which others selected clusters 
for X-ray observations, and those we imposed in selecting 
targets for GO observations. 

Our original intracluster light catalog comes from two 
sources and so itself has varied selection criteria: Gon- 
zalez et al. (2005) chose their clusters to have a single 
dominant brightest cluster galaxy (BCG) with a major 
axis position angle within 45° of cast-west (to accommo- 
date their drift-scan exposures) and no visual evidence 
of an ongoing merger: Kelson et al. (2008) chose their 
clusters to have a BCG that is kinematically stationary 
relative to the cluster, to appear relaxed in X-rays (when 
such data were available), and to be at sufficiently large 
redshift to not overfill the spectrograph slit length and 
close enough so that the kinematics of the low surface 
brightness ICS could be measured. 

Although we are unable to reconstruct the selection 
criteria for the X-ray archival sources, there is a natural 
tendency in the archival data to sample more X-ray lumi- 
nous systems and Bautz-Morgan (BM) (Bautz & Morgan 
1970) Type I clusters. For our GO observations, our aim 
was to observe more of the Gonzalez et al. (2005) clus- 
ters, for which clear distinctions between the BCG and 
ICS components has been measured, allowing us to make 
better optical/X-ray comparisons. 



The selection criteria described above do not lead to 
a representative sample of clusters. Instead, all clusters 
have a central BCG (BM type I or I-II with Abell 3705 
being the sole BM type III), Abell cluster richness class 
values of 1 or 2 (Abell 1958), and show no clear opti- 
cal evidence of an ongoing major merger. The clusters 
are within the redshift range 0.02 < z < 0.13 and have 
line-of-sight velocity dispersions, cr, spanning 500 to 1000 
km s~^, corresponding to the masses of rich groups and 
clusters (< lO" - 10^^ M©). We determine the size of 
these clusters by associating r2oo, the radius within which 
the mass overdensity is 200 times the universal value, 
with the virial radius (Cole & Lacey 1996). In prac- 
tice, we measure r2oo only for the Gonzalez et al. (2005) 
clusters, for which we have measured galaxy overdensi- 
ties (see Gonzalez et al. 2007, for details). We then use 
the theoretical relation between r2oo and a (Finn et al. 
2005), renormalized to match the Gonzalez et al. (2007) 
r200 values, for the Kelson et al. (2008) clusters. Wc use 
the rsoo — cr relation^ derived in Gonzalez et al. (2007) to 
compute the rsoo values for all our clusters. We present 
the physical parameters of the clusters in Table 1. The 
quoted cluster temperatures are the maximum annularly- 
averaged temperatures observed in our spectral fits (see 
§4.6 for details). This choice mitigates the ambiguity of 
the global temperature for cool core clusters and better 
reflects the depth of the cluster potential. 

3. DATA AND ANALYSIS 
3.1. X-ray Data: Preprocessing 

For our sample of 12 clusters, we have XMM-Newton 
EPIC X-ray (0.1 — 12 keV) observations, optical photom- 
etry (?'-band from Gonzalez et al. (2005) or R-band from 
Kelson et al. (2008)), and optical spectroscopy of mem- 
ber galaxies (Zaritsky et al. 2006). In this section, we 
focus on the analysis of the X-ray data. 

We present in Table 1 the observation ID of each X-ray 
data file, the date of observation, the source of the data, 
the filter used, EPIC-MOS camera exposure times, and 
the radius, router, out to which we perform spatially re- 
solved spectroscopy. The exposure times quoted are the 
average exposure times of the MOSl and M0S2 imaging 
spectrographs. We provide both the original imfiltered 
exposure time tMOS and the final fiare-filtered exposure 
time tMosji for each exposure. All MOS exposures were 
taken in Full Frame mode. 

We only analyze the MOS data because XMM-ESAS, a 
quiescent background modeling software that effectively 
models the time-variable instrument and particle back- 
ground spectra, is only available for MOS and because 
the XMM reflection grating spectrograph (RGS) does not 
offer the necessary fleld-of-view or sensitivity. We use 
XMM Science Analysis Software (SAS) 6.5 and FTOOLS 
5.3.1 to reduce the MOS data and create any XMM cal- 
ibration or ancillary data. 

The MOS data are prepared for further analysis as 
follows. We run the EMCHAIN task on the original data 
products to generate the event lists for each observation. 
We use recent calibration files (January 2006), which take 

^ All clusters in this paper lie in the a range (~ 500 to 1100 km 
s^-'-) where the rsoo ^ relation is calibrated. We advise against 
extrapolating this relation to lower velocity dispersions (cf. Balogh 
et al. 2008) 
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into account both the spatial and temporal variability of 
the MOS response (Stuhlinger et al. 2006). We choose all 
events, without any restriction on the energy range, with 
(PATTERN < 12) to filter out cosmic rays and other 
non-X-ray events. We set ((FLAG & 0x766a0f63) == 
0) to remove any bad pixels and retain the events that 
fall outside the ficld-of-vicw, namely in the unilluminated 
"corner" pixels used later for the background analysis. 

3.2. Background Subtraction 

Various sources contribute background to the observa- 
tions: 1) soft proton flares, which render the data co- 
incident with the flare useless, 2) a quiescent low-energy 
particle background, which exists even when no flares are 
present and is independent of the observed field, and 3) 
a "cosmological" background, which consists of resolved 
and unresolved sources in the field of view (for a complete 
review see Snowden et al. (2008)). 

XMM-Newton observations are plagued by fiaring 
problems from the satellite's passage through the Earth's 
radiation belt and proton clouds. These flares are typ- 
ically caused by soft protons that are tunneled into the 
telescope through the focusing optics and impinge on the 
detector. We remove these transient events before fur- 
ther analysis by excising all events for which the 2.5 to 
8.5 keV count rate exceeds a specified threshold within a 
time interval of 1 second. We determine the appropriate 
threshold value by generating a histogram of the count 
rate, fitting a Gaussian to the distribution to obtain the 
mean and dispersion adr, and identifying the count rate 
that is 2.5iTctr above the mean. We remove events in 
time intervals where the rate exceeds this threshold and 
iterate this procedure twice. The final integration times 
after fiare-filtering are given in Table 1. This filtration 
is fairly rudimentary because it is insensitive to low in- 
tensity proton flares. We account for this transient com- 
ponent in our spectral fitting (see §3.4 for details). In 
general, we find > 50% of the integration time of any 
given X-ray observation is useable for further analysis. 

The second background source consists of X-rays gen- 
erated within the instrument by the constant flux of low 
energy protons impinging on the detector and the sur- 
rounding housing. This background, both in the form of 
spectral lines due to fluorescence (Al and Si Ka) and 
continuum, changes temporally and spectrally, and is 
present even when no flaring is observed. All observa- 
tions must be corrected for it. We use the techniques out- 
lined by Kuntz & Snowden (2008) and implemented in 
the XMM-ESAS package for the quiescent particle back- 
ground subtraction. 

The procedure outlined by Kuntz & Snowden (2008) is 
a step-up in sophistication from typical quiescent back- 
ground subtraction techniques, which simply scale filter- 
wheel closed (FWC) data to match the corner pixel (i.e., 
dark) count rates. The quiescent particle background 
is known to vary temporally in spectral shape and in 
strength over long time scales (i.e., over several orbits). 
Corner pixels, which arc not exposed to photons because 
they lie outside of the field-of-view, are c;xposed only to 
this background over the same time interval as the ob- 
servation. Kuntz & Snowden (2008) characterize the ob- 
served corner pixel spectral using two spectral parame- 
ters: the high energy (2.5 to 12.0 keV) power law slope 
and the hardness ratio, which is the 2.5 — 5.0 keV to 



0.4 - 0.8 keV band flux ratio. 

There arc too few corner pixel counts to generate a 
background spectrum with sufficient S/N. Their method 
uses a database of FWC observations and extracts corner 
pixel data with similar slopes and hardness ratios for each 
MOS chip to augment and improve the S/N of the cor- 
ner pixel spectrum of the observation. The background 
spectrum within the region of interest is created by mul- 
tiplying the corner spectrum derived for each chip by the 
ratio of FWC object region spectrum and the FWC cor- 
ner pixel spectrum. Lastly, for a specified source region 
they combine the individual background components de- 
rived for each of the 7 chips using appropriate weights 
based on each chip's coverage of the source region and 
create a single background spectrum. 

We discuss our methods for addressing the "cosmolog- 
ical" background in the next section. 

3.3. X-ray Data: Spectroscopy 

To measure abundance and temperature profiles, we 
generate spatially resolved spectra for circular annuli 

centered on the clusters. Wc first identify the cluster 
center, set the annuli, locate point sources within each 
annulus, and remove point sources. We find the cluster 
center on a surface brightness map. Then, we radially 
bin the data in multiple annuli of variable widths rang- 
ing from 30 to 300 arcsec to obtain the highest radial 
resolution while maximizing the signal-to-noise in each 
bin. We choose 30 arcsec as a lower limit to minimize 
PSF blending effects between adjacent bins. The bin- 
ning choice varies for each cluster and is a function of the 
signal-to-noise of a particular observation. To remove 
point sources, we generate an image of each exposure 
and use the CIAO (Chandra Interactive Analysis of Ob- 
servations) 3.3 wavelet source detection utility, WAVDE- 
TECT. For each source, we obtain from WAVDETECT 
the elliptical aperture that contains 99.73% (within 3cr 
of a Gaussian distribution) of the source light. Elliptical 
apertures are necessary due to the off-axis changes in the 
XMM-Newton PSF. We visually inspect each source list 
and associated elliptical apertures to remove any spu- 
rious detections and ensure that the source extraction 
regions are large enough. We generate a region file con- 
taining the acceptable (point-source free) pixels in each 
bin. 

Next, we extract spectra and correct for the particle 
background spectrum. Using scripts available in XMM- 
ESAS that implement in the Kuntz & Snowden (2008) al- 
gorithm for determining the particle background, we gen- 
erate particle-background spectra that we subtract from 
our annularly binned cluster spectra. We compare the 
background and observed spectrum within an annulus to 
ascertain the quality of fiare-filtering. At high energies 
(>10 keV) where the telescope sensitivity is effectively 
zero, we do not in any case find that the background 
spectrum significantly deviates from the observed spec- 
trum, indicating that the fiares were effectively removed. 
If there remains low-level flaring, which cannot be de- 
tected by the above technique, we remove it by explicitly 
modeling it in our spectral fits as discussed below. 

Using GRPPHA, we group all of the response and 
background files together per spectrum, and we bin both 
the observed and background spectra to have 50 counts 
per energy bin, which ensures Gaussian statistics. Ide- 
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ally, we would want to construct energy bins to maximize 
signal-to-noise, but GRPPHA docs not have this feature. 
Count-binning is adequate for bins that are source count 
dominated, but inadequate for annuli that are particle 
background dominated. If the background counts are 
>25% of the signal counts, the high energy tail of the 
background-subtracted spectrum has very few counts, 
which results in a noisy spectrum. 

We use XSPEC 11 to fit physical models to the data 
and obtain measurements of temperature and chemical 
abundance. We fit to the 0.5 — 10 keV region of the 
spectrum, because it covers the Fc L-shcll 1 kcV) and 
Ka (6.7 keV) transitions, which are useful in determining 
the Fe abundance. XMM-Newton is known to have cali- 
bration issues at lower energies (Kirsch et al. 2005) and 
there are too few counts at energies greater than 10 keV. 
This range is required for proper determination of the 
continuum level and for characterizing the instrument 
background (Snowden et al. 2008). We use a complete 
model that accounts for the source, sky background, and 
instrumental emission, and fit all of the annuli simultane- 
ously. Our model of the sky background includes emis- 
sion from the local hot bubble, a residual soft proton 
flare background, and unresolved X-ray sources. We as- 
sume the cosmic background is constant over the field-of- 
view. Although this assumption is not entirely correct — 
the galactic X-ray background is patchy over scales of a 
few arcmin — it is reasonable because the XMM-Newton 
fleld-of-view is 30 arcmin, which effectively averages out 
small-scale variations. To constrain the low energy end of 
the background, we use additional ROSAT All Sky Sur- 
vey (RASS) data, which covers the 0.1 — 2.4 keV spectral 
range. We use the X-ray background tool at HEASARC 
that culls the RASS database to obtain a ROSAT posi- 
tion sensitive proportional counter (PSPC) spectrum of 
the annular region 1 to 2 degrees from the center of the 
cluster and the associated PSPC response file. Most of 
the emission in the 0.1 — 2.4 keV range is "background" 
because the cooler galactic plasma (0.1 keV) emits more 
efficiently at these energies than the cluster plasma (^ 
2 to 10 keV), therefore the PSPC spectrum provides a 
critical constraint to the soft cosmic background. 

We now fit the multicomponent model and follow to a 
large extent the prescription outlined in Snowden et al. 
(2008). The actual fitting procedure begins with the ap- 
plication of the IDL findspeclim routine to find the upper 
and lower spectral channels that correspond to the up- 
per and lower energy limits (0.5 — 10 keV). We then run 
genfitfile.pl, with known column density (Nh) derived 
from the Dickey & Lockman (1990) HI survey (using the 
HEASARC nH tool) and the cluster redshift from the 
NASA Extragalactic Database (NED) or our own data 
(Gonzalez et al. 2007), to generate an XSPEC initial 
fit file. We simultaneously fit both the MOSl and M0S2 
data as well as the PSPC data with XSPEC using the fol- 
lowing XSPEC model expression: "model bknpower/b + 
gaussian + gaussian + constant( apec + ( apec + apec + 
powerlaw )wabs + ( apec )wabs". bknpower/b \s a, hvoken 
power law, which is independent of the telescope's col- 
lecting area, that models any residual soft proton fiares 
and takes the following functional form: 



where f{E) represents the flux of the model as a func- 
tion of energy E, n is the normalization constant, a and 
P are power law exponents, and Eq is the break energy. 
We flx the break energy at 3.3 keV and allow the ex- 
ponents on either side of the break to vary freely along 
with the model's flux normalization. We assume that 
the soft proton flux is constant over the instrument, and 
we use the data from all annuli scaled appropriately by 
the annular area (in units of arcmin^) to constrain these 
parameters. The two Gaussians, signifled by gaussian 
in the XSPEC model expression above, represented in 
the analytical expression below model the instrumental 
X-ray florescence lines: Al Ka (1.49 keV), which is espe- 
cially apparent in annuli with few cluster counts, and Si 
Ka (1.74 keV). The model flux is then 



fiE): 
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E — En 
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f{E) = KE-°' E<Eq 
= kE^-''E-^ E>Eo 



(1) 



where k is the normalization constant, Eq is the line cen- 
ter, and a is the line width. Kuntz & Snowden (2008) 
have shown that the strength of these lines vary with 
each chip in the MOS instrument, so we flt for the flux 
normalization of these lines for each annulus while keep- 
ing the line energies and widths fixed. We initially also 
allowed the widths to vary but the best fit width val- 
ues were < 0.01 eV, which is less than the MOS spectral 
resolution. 

The remainder of the parameters model the cosmic 
background and the source, apec (Smith et al. 2001) 
is a collisional plasma model with the following free- 
parameters: temperatm'e (T), chemical abundance (Z), 
redshift (z), and flux normalization (k). wabs (Morri- 
son & McCammon 1983) implements a neutral hydrogen 
column photo-electric absorption model and is parame- 
terized by Nh- The constant accounts for the solid angle 
each annulus subtends on the sky (in arcmin^). For each 
annulus, we flx the value of the solid angle by appropri- 
ately scaling the value computed by SAS BACKSCALE, 
a utility that counts the number of square pixels within a 
given detector area. This is necessary to equally weight 
the fltting of background parameters over all annuli. The 
background is modeled with 4 components: the local 
unabsorbed hot bubble (apec, T ^ 0.1 keV, Z = Zq), 
the cool absorbed galactic halo {apec, T = 0.1 keV, 
Z = Zq, usually a negligible component), the warm ab- 
sorbed galactic component {apec, T ^ 0.3 keV, Z = Zq), 
and an unresolved extragalactic background. The extra- 
galactic background is modeled as a power law in en- 
ergy with an index of 1.41 and a flux normalization of 
8.88 X 10"'^ counts sec~^ keV"-^ arcmin"^. We fit the neu- 
tral hydrogen column of the absorbed models using the 
HI column density observed at the location of the cluster 
(Dickey & Lockman 1990) as an initial guess. Lastly, the 
cluster emission itself is modeled with an apec model. For 
each annulus, all parameters for both MOSl and M0S2 
spectra are tied together. We fix the redshift of the clus- 
ter to that in NED or to our own measured value (Gon- 
zalez et al. 2007). We adopt Anders & Grcvcsse (1989) 
solar photospheric abundances. We use minimization 
to obtain best fit values. The fits provide a measure of 
the cluster temperature, Fe abundance, and flux normal- 
ization for each annulus. The maximum radius of our 
spectral fits, Vspec, given in Table 1 is determined to be 
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the radius within which the fractional error of the derived 
cluster abundance does not exceed 50%. 

3.4. X-ray Data: Imaging 
3.4.1. Generation of Smoothed Maps 

Our principal objective for the imaging is to exam- 
ine the cluster morphologies, generate surface brightness 
profiles, and perform surface brightness fits. Wc take the 
MOSl and M0S2 flare-filtered event lists of each target 
and generate a soft (0.3 — 1.25 keV) and a hard (2 — 8 
keV) image. We choose these energy ranges to avoid 
contamination from the Al and Si Ka fluorescence lines. 
We generate exposure maps with EEXPMAP for each 
image, which are used to flatten the data and remove vi- 
gnetting. We create bad pixel masks using EMASK. We 
identify bad columns by searching for steep gradients in 
the exposure map and remove them using the gradient 
parameter in EMASK. Using the source region files pro- 
duced for spectral fitting, wc mask out point sources. 
For each image, we use XMM-BACK-IMAGE to gen- 
erate a corresponding particle background image. We 
mosaic the background-subtracted images to form a sin- 
gle 0.3 — 8 kcV image. If the data at a given sky position 
are available in multiple exposures, we add the fluxes. If 
data are available at a given sky position from only one 
pixel, we scale the flux up to the total exposure time. If 
data are not available at a given sky position, i.e., if the 
corresponding pixel is masked in both MOS exposures, 
we set the image value to zero. 

To produce higher signal-to-noise images, we adap- 
tively smooth the mosaicked image using ASMOOTH. 
Our goal is signal-to-noise ~ 20 over the cluster, and so 
wc choose a maximum kernel width of 10 pixels (25 arc- 
sec), no minimum size, and ignore image values of zero. 
We present the smoothed images in Figure 1. Contours 
are spaced logarithmically and the lowest contour is set 
at a fixed value of 2.0 x 10~^ counts s~^ arcsec"^ chosen 
to be unaffected by the noise present in the smoothed 
image. 

3.4.2. Surface Brightness Fitting 

To measure and model the Fe mass, we require a mea- 
surement of the gas density profile. We derive this den- 
sity profile from the 2D surface brightness profile fitting 
of MOS images and the spectral fits. The surface bright- 
ness profile constrains the shape of the gas density distri- 
bution, whereas the spectral fits provide the normaliza- 
tion for that distribution. We use SHERPA, a versatile 
fitting package that is part of the CIAO 3.3 distribu- 
tion, to fit the 2D model. Starting with the particle- 
background-subtractcd, exposure-corrected images, we 
prepare these images for fitting by first adding a small off- 
set of 1.6 X 10~^ counts s~^ arcsec"^ for the soft band and 
3.2 X 10~^ counts s~^ arcsec"^ for the hard band, which 
excludes negative pixels. We then multiply the images 
with an unvignetted exposure map and add the particle 
background to the image. For this fit, wc use the unvi- 
gnetted exposure maps with point source regions zeroed 
out, and convolve the model with the 2D point spread 
function (PSF) generated by CALVIEW. Because we are 
fitting an extended source, off-axis PSF effects are not 
as important, so we use the on-axis PSF over the en- 
tire field. For simplicity and consistency with previous 



studies, we choose to fit the canonical beta model (Cava- 
liere & Fusco-Femiano 1978), which is a straightforward 
relation between surface brightness Ix{R) and volume 
density Pg{r): 



and 
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where the free parameters /x(0) and Pg(0) are model 
normalizations, rx is the core radius, and /3 defines the 
rate at which these profiles decline with projected radius 
i? and unprojected radius r, respectively. 

We model the background, which includes both the 
vignetting-corrected sky and unvignetted particle back- 
ground, with a constant term. This approach produces 
consistent results with fits performed using blank-sky 
subtracted images in two representative fields. For the 
blank-sky subtraction, we use the Carter & Read (2007) 
blank-sky files that were observed though the same filter 
as our observations. These files include both the parti- 
cle and cosmic backgrounds. We cast the blank-sky event 
list files into our sky coordinates and scale them to match 
the exposure time of our exposures before carrying out 
the background subtraction. 

For our fits, we simultaneously fit both the soft and 
hard band data, using the respective unvignetted expo- 
sure maps. Wc use the Cash statistic, which is appropri- 
ate at low count levels, to evaluate goodness-of-fit and 
the simplex algorithm to find the best fit parameters. 
We look at residuals to determine if there are clusters 
with bright cores. In clusters with unusually bright cen- 
ters (sec Appendix for more details), we mask the cen- 
tral ~ 30" and refit the data. The best fit for Abell 
496 converges to one with zero background flux, which is 
unphysical. We attribute this to large residual non-axis 
symmetric features found in the inner 4' of the cluster 
(see Figure 4). Such anomalies are rare and are discussed 
in the Appendix. 

Initially, using a simple spherically symmetric beta 
model, we fit for the center, core radius, and constant 
background level. The best fit had large residuals in 
most cluster centers, (see Figure 3 for one representative 
case), which suggested that we consider more complex 
models. Therefore, we included variable ellipticity and 
position angle. This change improved the fits, and the 
resulting residuals were no longer globally asymmetric; 
however, the statistical improvement was marginal. Be- 
cause cluster X-ray emission can have two temperature 
components (a lower temperature, higher density com- 
ponent near the core, and a higher temperature, lower 
density component at larger radii; Jia et al. (2004)), we 
then explored fitting two elliptical beta models. The two 
component model performed only marginally better than 
the single component model and had poorly constrained 
outer component parameters as well as /? values that were 
unphysically large (> 2). We therefore adopt our single 
component, spherical beta model fits to maintain consis- 
tency with our spectral analyses, which use circular and 
not elliptical annuli, and to simplify our mass calculation. 
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To account for the fact that wo arc observing surface 
densities but want to work with volume densities, we 
must deproject the observations to normalize the den- 
sity model. To do this deprojection, we model the cluster 
as spherically symmetric with a gas density profile that 
obeys Equation 4. We calculate the gas density using 
the Abell transform and compute the central gas den- 
sity that reproduces the derived cluster apec normaliza- 
tion constant k (given below) for each annulus, assuming 
hh = 0.82ne (Pizzolato et al. 2003): 

"^ 4.p^(i+z))^ y"-""^^' 

where Da is the angular diameter distance (in cm), z 
is the redshift, Ue and njj are the electron density and 

the Hydrogen density in cm~^, respectively. We average 
the central gas densities derived from fitting to individ- 
ual annuli for use as the central gas density Ph{0)- We 
compute masses by integrating Equation 4. 

3.5. Monte Carlo Error Modeling 

We determine all central density and cluster mass er- 
rors through Monte Carlo simulations. We adopt the 
spectral and radial profile best fit values and the la errors 
derived by XSPEC and SHERPA. The error distribution 
is modeled as Gaussian. SHERPA^s projection function 
carries out its own proper Monte Carlo simulation for a 
given fit and returns asymmetric la errors. We model 
the asymmetric errors with two Gaussians. We carry out 
2000 trials per cluster where for each trial we generate 
new values of spectral normalization, abundance, f3, and 
rx, and compute the masses. The resulting distribution 
of derived masses is slightly asymmetric, but centrally 
peaked. Increasing the number of trials does not change 
the resulting distribution. Using the resulting distribu- 
tion of mass estimates, we evaluate our la confidence 
limits. 

3.6. X-ray Luminosities 

Wc measure cluster X-ray luminosities by combining 
the beta profile and spectral fits. Therefore, we can only 
do this for clusters where we were able to fit a beta pro- 
file, which excludes Abell 2877 and 3705. We first com- 
pute the luminosity enclosed within the radius of our 
spectral fit by using the normalization parameters of the 
apec model for each annulus. To determine the luminos- 
ity lying outside of the spectral fit radius, we extrapo- 
late the luminosity out to r2oo using the beta model fit. 
To accurately quantify the errors in this luminosity, we 
carry out a Monte Carlo simulation using XSPEC 12. 
We first load the best fit spectral model into XSPEC 
and compute new parameters drawn from a multivariate 
gaussian for each trial. We then use these new param- 
eters to derive the new luminosity of the cluster within 
the spectral fit radius after zeroing out all background 
components. This procedure adequately accounts for the 
error in the complete spectral model rather than in just 
a simple model considering only cluster emission. We 
also include the errors in our surface brightness fit in our 
Monte Carlo simulation when we compute the total X- 
ray luminosity within r2oo- We carry out a total of 200 
trials and tabulate our results in Table 1. 



4. RESULTS 
4.1. Cluster Morphology and Profiles 

The X-ray cluster emission from our sample clusters 
is generally smooth, mildly elliptical, and centered on 
the BCG. In Figure 1, we show the 0.3 - 8.0 keV X- 
ray contours overlaid over the corresponding Digital Sky 
Survey"^ (DSS) field. AbcU 3693 has a clump to the 
southeast, although the central cluster appears relaxed. 
Abell 2877 and 3705 are exceptions that do not appear 
to be dynamically relaxed. Abell 2877's diffuse emis- 
sion is not centered on the BCG. Moreover, Abell 2877's 
low signal-to-noise and a bright, nearly central (2 ar- 
cmin from the cluster center) point source complicates 
the analysis. The NED database indicates that A2877 
consists of two superposed clumps, which further sug- 
gests that it is unrelaxed. Abell 3705 has asymmetric, 
unsmooth contours and a clump to the northwest, and it 
is our sole BM Type III cluster. Because we model our 
clusters as relaxed and symmetric objects, the results for 
these two clusters are suspect. 

In the left panel of Figure 2, we plot the ID surface 
brightness profiles with the point sources and particle 
background removed. Error bars represent the standard 
deviation of the mean surface brightness value for a given 
bin. Our results from fitting the single 2D beta model fits 
are in Table 1, and the radially averaged surface bright- 
ness fit residual profile (in units of x) is presented in a 
sub-panel beneath each cluster's surface brightness pro- 
file. We also compute ID values using our best fit 
2D model, which serve as an additional measure of the 
goodness-of-fit. Our fits to Abell 2877 and 3705 yield 
unphysical (3 and rx- We do not quote errors for Abell 
2877 fit parameters because SHERPA's Monte Carlo er- 
ror estimator fails to converge on a result. Because of the 
evidence suggesting that Abell 2877 and 3705 are unre- 
laxed, we do not include either in our final BCG-I-ICS 
versus cluster galaxy enrichment comparisons. 

4.2. Cluster Spectral Fits 

We plot the results of our spatially-resolved spectral fit- 
ting in Figure 2 with temperature profiles in the middle 
panel and chemical abundance profiles in the right panel. 
Error bars denote la values derived using XSPEC. Sev- 
eral of our clusters have cool cores and steeply declining 
abundance gradients that peak at ~ T^Fe,Q and fall to 
the canonical value of 0.3 T^Fe,® at large radii. These re- 
sults illustrate the importance of spatially-resolved spec- 
troscopy out to large radii. Only four of our clusters, 
Abell 1651, Abell 2811, AbcU 3693, and Abell 3705, do 
not show steep abundance gradients. Although the con- 
ventional wisdom is that clusters with isothermal profiles 
typically do not have abundance gradients (De Grandi & 
Molendi 2001), we see here that even clusters that appear 
isothermal, AbcU 2984, S84, and S540, have abimdancc 
gradients. In addition to abundance gradients, a few of 
our clusters show declining temperatures at large radii. 

® The Digitized Sky Surveys were produced at the Space Tele- 
scope Science Institute under U.S. Government grant NAG W- 
2166. The images of these surveys are based on photographic data 
obtained using the Oschin Schmidt Telescope on Palomar Moun- 
tain and the UK Schmidt Telescope. The plates were processed 
into the present compressed digital form with the permission of 
these institutions. 
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This may be an artifact of the spectral fitting arising 
from an over-subtracted background, but similar declines 
are observed by Vikhlinin et al. (2005) in their Chandra 
study of multiple clusters using a different analysis, and 
by Snowden et al. (2008), which suggests a physical ef- 
fect. 

4.3. Intracluster Gas and Fe Masses 

We use the following prescription to calculate the de- 
projected gas mass within any specified radius, which we 
determine by integrating the volume gas density profile 
(Eq. 4). We choose to evaluate quantities to a radius of 
r5oo, the radius at which the cluster mass density is 500 
times the Universe's current mass density and the largest 
radius for which gas and stellar mass measurements are 
available for most clusters. Reaching this radius often 
requires some extrapolation, so we also quote gas mass 
values out to the maximum radius of our spectral fits, 

^spec • 

Using the volume gas density profile, we compute the 
unprojected hydrogen gas mass enclosed within each 
annulus used in our spectral fits, multiply it by the 
spectrally-derived abundance to obtain Fe mass, and 
then use the solar photospheric Fe/H mass fraction 
(0.0026; Anders & Grevcsse 1989) to determine the phys- 
ical quantity of Fe in solar masses present in that radial 
bin. For r > ispec, we assume a constant abundance 
equal to that at Vspec- The recently revised solar values 
(Asplund et al. 2005) do not affect our results because 
our spectral fits also use the Anders & Grevesse (1989) 
solar photospheric abundances and the derived Fe abun- 
dances should thus scale linearly with the difference be- 
tween the two solar abundance catalogues. We present 
the calculated hydrogen and Fe masses in Tables 1, 1, 
and 1. Typically, the extrapolation is less of an issue 
for gas masses than for Fe masses because the siirfacc 
brightness fits that constrain the shape of the gas den- 
sity extend farther than igpec- Gas masses provided in 
these Tables are for hydrogen (H) because abundances 
are quoted relative to H by convention. To obtain the 
total gas mass, one must accoimt for the mass fraction 
of helimn and midtiply the hydrogen gas mass by 1.33 
(assuming primordial relative abundance). This adjust- 
ment is done in §4.5 for the comparison of our gas mass 
estimates with those presented in other studies. 

4.4. Systematic Error arising from Spherical 
Symmetry Assumption 

One possible source of systematic errors that could af- 
fect both our and previous studies is the assumption of 
spherical symmetry. Even relaxed looking clusters are 
often elliptical (Figure 1). Because the gas mass predic- 
tion is critically linked to the assumed geometry of the 
system, we test how the spherical assumption affects our 
results by carrying out the analysis on gas distributed 
in different ellipsoidal geometries. We start with the 
model gas density prescription outlined with (3 = 2/3 
and rx = 5 in §3.4.2, and recast it into an elliptical 
form: 

-3/3/2 

y 



Pg{x,y,z) oc 



1 + 



' X 



+ 



+ 



(6) 



where the e's arc the ellipticities along each of the three 
principal axes. We calculate the gas density on a 401'^ 
grid with the line-of-sight along the z-axis. We gener- 
ate five different models: a spherical one, two oblate 
spheroids with ellipticities (e^; = l,ey = 0.7, €z = 1) 
and (1,1,0.7), respectively, and two prolate spheroids 
with ellipticities (1, 0.7, 0.7) and (0.7, 0.7, 1), respectively. 
We also bound all of the distributions such that {{x — 

Xcen)/(^xf + ((y - yceu)/<iyf + {{z - Zcen) / < r^ax 

where Vmax = 200. Because the X-ray emission is thermal 
bremsstrahlung, we compute a surface brightness map by 
assuming isothermality, squaring the density profile, and 
integrating along the line of sight. We generate radial 
profiles by azimuthally averaging the surface brightness 
map. 

Figure 5 shows that the spherical model has the high- 
est siirface brightness, followed by the oblate and then 
prolate cases. By normalizing the profiles at a given ra- 
dius, we see similarities in the projected profiles despite 
the different geometries; for example, all have similar 
slopes after the knee in the curve, which means that they 
have similar (3 values. Starting with a cluster profile as 
shown in the right panel of Figure 5, we compute the gas 
mass for each ellipsoidal geometry and compare to the 
spherical case. We determine the upper and lower error 
bounds for the oblate and prolate cases. For the oblate 
case, the error ranges from +19.5% to +11.5%, relative to 
the spherical case, where the largest error corresponds to 
the cluster being flattened along the line of sight. For the 
prolate case, the upper and lower bounds are +17.5% and 
+25.9%, where the worse case error arises when the clus- 
ter is flattened along the line-of-sight. We conclude that, 
except for extreme geometries, the error arising from our 
spherical assumption is likely to be 15%. 

4.5. Comparison with Previous Work 

As a check on onv procedures, we compare our temper- 
ature, luminosity and gas mass measurements to those 
from previous studies. We first compare our luminosity 
measurements with the REFLEX (Bohringer et al. 2004) 
values for the systems that we have in common. To do 
so, we recompute our 0.1 — 2.4 keV luminosity within 
the REFLEX aperture. In the left panel of Figure 6, we 
show the fractional difference between our respective lu- 
minosity measurements. In general, we agree to within 
20% and there is no systematic trend with cluster veloc- 
ity dispersion. There are a few significant outliers, such 
as Abell 2984 and 3693, for which our luminosities are 
lower by as much as 60%. There is no clear explanation 
for this discrepancy. We also compare the peak temper- 
atures for seven of our clusters with those found in the 
BAX database (Sadat et al. 2004) in the right panel of 
Figure 6. In general, the values are consistent. Our val- 
ues are systematically higher for four clusters (three of 
which have cool-cores), but this shift is likely to arise 
from our use of peak temperatures instead of emission- 
weighted temperatures, which are usually dominated by 
cool-corcs. The one exception is Abcll 2811, which has 
an isothermal temperature profile. 

We also compare our spectral fit values for the three 
clusters, Abell 496, 3112, and 4059, that are in the recent 
cluster catalogue released by Snowden et al. (2008) with 
temperature, abundance, and surface brightness profiles. 
In all cases, our profiles are consistent with those derived 
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by Snowden ct al. (2008). We do see a small systematic 
offset in the temperature profile of Abell 3112, but the 
difference is less than 10%. In general, Snowden et al. 
(2008) measure temperatures to larger radii, but they 
measure abundance values to similar radii as we do. We 
do not attempt to extend our temperature fits to larger 
radii as that is not the focus of this paper. 

Lastly, wc carry out gas mass comparisons, which re- 
flect on our ability to ultimately determine Fe masses. 
We search the literature for gas mass measurements 
within rsoo for our clusters (Reiprich 2001; Castillo- 
Morales & Schindler 2003; Piffarctti et al. 2005) and 
tabulate the results of our comparisons in Table 1. For 
the comparison, we determine the gas mass interior to 
the rsoo defined in each comparison study and convert 
to their cosmologies. We hereafter refer to Einstein- 
de Sitter cosmology {il„i = 1, ^^a = 0, Hq = 50 km 
Mpc^^) as SCDM50. The methods used in papers 
with ROSAT data (Reiprich 2001; Castillo-Morales & 
Schindler 2003) involve fitting a projected gas mass pro- 
file to the observed surface brightness distribution, with 
slight variations in the choice of fitting function, although 
use of the beta profile is common. Our gas masses are 
consistent to within our uncertainties for one of three 
clusters, Abell 3112, from Castillo-Morales & Schindler 
(2003). For Abell 496 and 1651, the results differ by 18% 
and 11%, respectively, but Castillo-Morales & Schindler 
(2003) do not provide error estimates for their gas masses 
so we cannot judge the significance of this disagreement. 
Our results agree with those of Reiprich (2001) for five 
clusters. 

We also compare our gas mass values with those of 
Piffaretti et al. (2005), who, like us, use XMM-Newton 
observations. Our masses are systematically higher than 
theirs by as much as 50% for the same LCDM70 cosmol- 
ogy. They employ a more sophisticated spectral depro- 
jection technique to derive their masses, yet our appli- 
cation of a deprojection technique to one of our clusters 
suggests that this makes only a modest change and is not 
responsible for the difference. Instead, the discrepancy 
arises from the gas density profiles. Our gas densities, in- 
dependent of radius, are consistently higher than theirs 
by a factor of ~ 1.6. 

If the Piffaretti et al. (2005) results arc closer to the 
correct value, which would suggest that all previous 
ROSAT measurements are incorrect, we are overpredict- 
ing the cluster Fe mass by 40 — 50% and hence over- 
stating any discrepancy between the predicted Fe mass 
created by the ICS, BCG, and cluster galaxies and the 
measured Fe mass in the ICM (§5.1). Therefore, we carry 
out an additional test by plotting the relation between 
seven other clusters that have gas masses measured by 
Vikhlinin et al. (2006) and velocity dispersions from the 
literature (see Gonzalez et al. 2007). We then compare 
this relation to ours (Figure 7) . We do not see a system- 
atic offset between the mean trends found in our data 
and the Vikhlinin et al. (2006) data, suggesting further 
that our gas masses are accurate. 

4.6. X-ray Scaling Relations 

As a further test of our spectral fitting techniques and 
the accuracy of the velocity dispersions of our clusters, 
we compare our sample to known Lx — T — ct— optical 
relations. We define the cluster temperature, which is a 



reflection of its dynamical mass, to be the highest tem- 
perature derived in our spatially- resolved spectral fits. 
This choice is reasonable because the derived tempera- 
tures are annularly-averaged and not likely to be affected 
by localized hot spots. In Figure 8, we show that the ma- 
jority of our clusters fall within the errors of the T — a 
relation given by Wu et al. (1999): 

^ _ -|^q2.49±0.032i0.64±0.02 j-j-j 

which is derived from a fit to 92 clusters. Outliers from 
this, and the other scaling relations discussed below, are 
described individually in the Appendix. 

As a further sanity check, we determine the Lx — ct 
and Lx — T relations for the 10 of our clusters with 
measured luminosities. We can adequately fit a power 
law in both cases (Figure 8), obtaining the following two 
relationships: 

Lx = 1036-ll±2.1V2-83±0.70 (g) 

and 

j^^ _ -[^Q42.38±0.46y2.79±0.77 j-g-j 

where Lx is the luminosity measured in the 0.5 — 2.0 keV 
range. These relationships establish the validity of the 
X-ray reduction and analysis techniques employed here. 

5. DISCUSSION 

There have been numerous previous attempts to ex- 
plain the metals in the intrachistcr medium (ICM) (see 
Portinari et al. 2004, and references therein). If cluster 
galaxies are the sole polluters of the ICM, which is the 
basic supposition of these studies, then it is not possible 
to produce the observed Fe mass using standard assump- 
tions for the IMF, SN la rate as a function of time, and 
metal loss efficiency of galaxies (Portinari et al. 2004). 
This has led to arguments for a top-heavy or a variable 
IMF (Matteucci & Gibson 1995; Gibson & Matteucci 
1997; Moretti et al. 2003; Portinari et al. 2004; Loewen- 
stein 2006), which would reduce the metal mass locked 
in low-mass stars with lifetimes longer than a Hubble 
time (the locked-up fraction) and produce a larger num- 
ber of Type II supernovae. These studies often assume 
that strong winds from the supernovae blow a significant 
fraction of the metals out of the galaxies, i.e., up to a 75% 
metal loss fraction from elliptical galaxies, which leaves 
behind the 25% locked up in lower mass stars. Recent 
work produces a more consistent picture of the ICM en- 
richment through hydro dynamical simulations that rely 
on galactic superwinds (Romeo et al. 2006) and possibly 
ram-pressure stripping (Kapferer et al. 2007) in order to 
explain the observed trend in Fe abundance with redshift 
(Calura et al. 2007). Although these various models pro- 
duce consistent results, the underlying prescriptions have 
yet to be observationally proven and/or calibrated. 

Often neglected in these analyses is the contribution 
from intracluster stars that pollute the ICM in-situ. 
With the measurements of the cluster Fe mass and the 
luminosity of ICS in hand, we are in position to address 
a variety of questions regarding the role of intracluster 
stars in the chemical enrichment history of clusters. Do 
intracluster stars play a significant role? If so, how does 
their contribution compare to that of stars in galaxies? 
Can the two components together account for the metals 
observed in the ICM? Can we use these constraints to 
place limits on the metal mass loss from galaxies? Do we 
require non-standard assumptions to explain the metals? 
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5.1. Metals from Intracluster Stars 

We begin the calculation of the metal contribution 

from intracluster stars using the parameters for the in- 
tracluster light measured by Gonzalez et al. (2005) and 
Kelson et al. (2008). We then calculate the mass-to- light 
ratio, corresponding SN rate, and chemical deposition 
rate using the PEGASE.2 evolutionary models (Fioc & 
Rocca-Volmerange 1997) for 1 Mq initial parcel of gas. 
Scaling those results by the current mass of intracluster 
stars and integrating the SN rate over the lifetime of the 
stellar population, we calculate the Fe mass produced 
by the ICS. Because decomposition of the light from the 
brightest cluster galaxy and the intracluster stars is am- 
biguous, particularly in the Kelson et al. (2008) data, 
which do not extend as far in radius as the Gonzalez et 
al. (2005) data, we compute the Fe mass generated by 
the combination of the BCG and intracluster stars and 
refer to that as the BCG+ICS contribution. Where we 
can decompose these two components because of abrupt 
changes in position angle or ellipticity. the BCG gener- 
ally contributes only 10-20% of the light (Gonzalez et al. 
2005). Therefore, even though we present results for the 
BCG and intracluster stellar components together, they 
do primarily reflect the metal contribution from intra- 
cluster stars. 

We make two important simplifying assumptions. 
First, we consider the ICS to be an old population orig- 
inating from an initial, short-lived episode of star for- 
mation. Williams et al. (2007), in their study of intra- 
cluster stars in Virgo, conclude that the majority of the 
ICS are old (> 10 Gyr). Second, we assume that any 
metals produced by the ICS find their way into the ICM, 
regardless of the origin of the ICS. Recent simulations 
suggest that the ICS are predominantly stripped from 
parent galaxies through tidal interactions (Murante et 
al. 2007). There are also observational claims that these 
stars form in gas being stripped from infalling cluster spi- 
rals (Sun et al. 2007). All the metals produced by stars 
that eventually become ICS, either before, during, or af- 
ter the tidal interactions, are assumed to lie in the ICM. 
If, for some reason, the interactions that created the ICS 
were less efficient at removing metals from galaxies and 
if much of the metal enrichment due to these stars hap- 
pened prior to their becoming ICS, then we will over- 
predict the ICS Fe contribution. The second assumption 
can be tested observationally by measuring present-day 
intracluster and cluster galaxy SN la rates. There are 
indications that ~ 20% of cluster SN la explosions occur 
outside of galaxies (Gal- Yam et al. 2003; Domainko et 
al. 2004), though these results are encumbered by large 
uncertainties. An extensive cluster SN la survey is un- 
derway to improve the precision of this result (Sand et 
al. 2008). 

We model the ICS as a 14 Gyr old stellar system with 

an initial, singi(\ 1 Gyr long, constant star formation rate 
burst and a standard Miller-Scalo initial mass function 
(Miller & Scalo 1979). We choose default PEGASE pa- 
rameters: Model B yields for SN II ejecta (Woosley &: 
Weaver 1995), no gas infall, no galactic winds, no neb- 
ular emission, and no extinction. We use 250 Myr time 
steps for the calculation. Using shorter, 1 Myr, time 
steps changes the final values only by a few percent. We 
set the initial metallicity to 0.004 by mass fraction. This 



is a plausible upper limit (0.2 Zq), which is ~ 2x larger 
than that of the old populations in local galaxies (e.g., 
the SMC (Harris & Zaritsky 2004) and LMC (Harris & 
Zaritsky 2008)). Taking a more conservative initial value 
of 0.0001 leads to a modest decline in Fe of 5% (largely 
attributable to the change in mass-to- light (M/L) ra- 
tio), which is well within the other micertainties. We 
do not use the evolving metal abundances calculated by 
PEGASE in any way other than indirectly through the 
calculated M/L. Hereafter, we refer to this model as our 
baseline model. 

The true nature of SN la, which are the dominant 
source of Fe, and their progenitors remains a subject of 
debate (Maoz & Mannucci 2008). Currently, two possible 
pathways are theorized to produce Type la supernovae. 
In the single-degenerate (SD) case, a single white dwarf 
(WD) accretes matter from a close companion until it ex- 
plodes upon reaching the Chandrasekar limit (Whelan & 
Iben 1973). The double-degenerate (DD) model involves 
the merger of two binary WDs (Iben & Tutukov 1984) 
that leads to the explosion. We model the SN la rate as 
a function of time by using the rates calculated by PE- 
GASE, which uses a popular SD prescription (Greggio & 
Renzini 1983). This prescription has a free parameter, 
A, the close binary fraction, that is manually tuned to 
match the present day SN la rate. We choose A = 0.175 
to match the observed present-day, cluster early-type SN 
la rate (0.066lj];!]^^ SNuM, Mannucci et al. (2008)). The 
steep evolution in our model rate as a function of look- 
back time is also consistent with all other cluster SN la 
rate measurements carried out at different epochs {z <1; 
Gal- Yam et al. 2002; Sharon et al. 2007). The purpose 
of this modeling is to provide the SN rate over the life- 
time of the stellar population. However, the dependences 
of the SN rate on environment and look-back time are 
poorly known. Adding to the uncertainties surrounding 
the SN rates, Scannapieco & Bildsten (2005) recently 
presented evidence for two classes of la's, prompt and 
delayed, whose relative importance will naturally vary 
over time and with environment. We consider a crude 
empirical model incorporating these two SNe la classes 
at the end of this section. 

We adopt empirically determined Fe yields and de- 
fine yield as the mass of Fe ejected into the interstellar 
medium per SN explosion. We use a Type la Fe yield 
of 0.7 Mq (Maoz & Gal- Yam 2004), which is consistent 
with both the modeling of bolometric light curves of SNIa 
(Contardo et al. 2000) and theoretical SN la yield calcu- 
lations (Thielemann et al. 1986). Given the large scatter 
in Type II yields (0.0016 to 0.26 Mq; Hamuy (2003)), 
we adopt an intermediate mean yield of 0.05 Mq from 
Elmhamdi et al. (2003), who estimate the mass of ^^Ni 
(which eventually decays to Fe) produced in a Type II 
SN by matching calculations to measured light curves 
and Ha luminosities. We make the instantaneous re- 
cycling approximation by adopting a constant yield per 
Type II SN. 

To scale the PEGASE results to the corresponding stel- 
lar mass of the ICS, we deproject the optical luminosity 
profile using a 3-D Hernquist profile (Hernquist 1990), 
which is known to fit a deVaucouleurs profile in pro- 
jected space. This procedure enables us to obtain the 
total optical light within spherical shells defined by the 
corresponding X-ray spectral fit annuli. The optical lu- 
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minosity is converted to a stellar mass using an i-band or 
R-band mass-to-light ratio of 2.25 or 2.95, respectively. 
The enclosed mass is given as 

^W = ^*°*(7T^' (10) 

where Mtot is the total mass derived from the projected 
surface brightness profile and a is related to the scale 
length, i?e, through the relation i?g « 1.8153a. 

First, we use the PEGASE model to calculate the Fe 
mass produced by the BCG+ICS and compare those val- 
ues to the observed ICM Fe masses. We do not con- 
sider the locked-up Fe fraction in the ICS because re- 
cent metallicity measurements of the ICS in Virgo give 
low values (~ 0.1 on average; Williams et al. 2007). 
Consequently, only a negligible fraction of Fe (~ 2%) 
produced by the ICS is locked up in stars and virtually 
all of the metals produced by the ICS go into the ICM. 
Our results for volumes inside radii of rgpec and rsoo are 
presented in Tables 1 and 1, respectively. Column 2, in 
both Tables, contains the measured Fe mass within the 
specified radius, while Columns 3, 4, and 5 contain the 
predicted Fe mass contributed by the BCG only, the in- 
tracluster stars only, and the combination of these two 
components, respectively. We are not able to decompose 
the BCG and ICS in all systems. The final column lists 
the fraction of the observed Fe mass contributed by the 
BCG+intracluster stars for the baseline model. To esti- 
mate the uncertainties in the fractional contributions, we 
adopt a 15% systematic error in the BCG+ICS Fe con- 
tribution arising from variations in M/L and from pro- 
jection effects. We incorporate the random error from 
X-ray Fe mass measurements. We do not include the 
uncertainty from the SN la rate in the individual clus- 
ter values, but include it in the averages that we discuss 
below. We present the r^pec results for the sake of com- 
pleteness but do not use them in our work. 

We place all clusters on a common spatial scale by cal- 
culating all of the relevant quantities out to rsoo. Becaiise 
the ICS are more centrally distributed than the metals, 
this choice of radius is important. We chose r5oo because 
it requires the least amount of extrapolation from Vspec 
(see Table 1) and includes a significant portion of the 
cluster metals (cf. Ettori (2005), Calura et al. (2007)). 
All of our clusters' Fe profiles are either flat or decreas- 
ing with radius (Figure 2). Therefore, our method of ex- 
trapolation will at worst overestimate the amount of Fe 
present within rsoo in clusters and thereby lead us to un- 
derstate the contribution from intracluster stars within 
that radius. 

Our results indicate that intracluster stars do play a 
significant role in the enrichment history of clusters. In 
our baseline model, the BCG-I-ICS contribute an average 
of 31^9^% of the observed intracluster Fe within rsoo. We 
show in Figure 9 (top panel) the fractional contribution 
of the BCG-I-ICS for each cluster with the dotted line 
representing the 31% weighted average value. All aver- 
ages are computed using 1/error^ as weights. The highly 
discrepant point with the lowest velocity dispersion is 
Abell 2984 and incidentally has the largest errors in frac- 
tional BCG-I-ICS contribution. The rsoo values for two 
lowest velocity dispersion clusters (Abell 2984 and Abell 
S 84) are only constrained by a single low velocity dis- 
persion point in the calibration of the rsoo — f relation 



(Gonzalez et al. 2007), and we suggest some caution in 
interpreting the results from these two clusters. Because 
the BCG typically contributes much less Fe than the ICS 
in clusters where we can separate the BCG and ICS, we 
conclude that the ICS are significant polluters of the in- 
tracluster medium within rsoo and need to be included 
in all future ICM metal budget accounting. 

Next, we explore the potential effect of having two 
classes of SN la's. We restrict ourselves to the parame- 
ters of the two populations as specified by Scannapieco & 
Bildsten (2005). Their model consists of a prompt popu- 
lation (with explosion delay time r ^ 0.5 Gyr) that scales 
in number with the star formation rate and of a delayed 
population (r ^ few Gyr) that scales with the total stel- 
lar mass. More luminous (prompt) SNe la are observed 
preferentially in star- forming galaxies, while fainter (de- 
layed) ones are seen in red galaxies with little star forma- 
tion (Scannapieco & Bildsten 2005). Using the prescrip- 
tion suggested by Scannapieco & Bildsten (2005), we de- 
rive the SN la rate using the star formation rates and stel- 
lar masses output by PEGASE for a 14 Gyr population. 
Most of the Fe enrichment is dominated by the prompt 
component at early times when the ICS progenitors are 
actively forming stars. We obtain a BCG-I-ICS average 
fractional Fe contribution of 22lg% within rsoo, which 
is consistent with our one-component, baseline model. 
In contrast with the baseline model, the two-component 
model is relatively insensitive to the present day SN la 
rate (rescaling this delayed population upward by a fac- 
tor of 3 increases the fractional Fe contribution by < 5%, 
because the prompt population is the dominant Fe con- 
tributor). The Scannapieco & Bildsten (2005) prescrip- 
tion for the prompt rate is based on measurements of 
field galaxies over a range of redshifts. Because we do 
not know how the prompt component scales with envi- 
ronment, we retain their scaling. Although the baseline 
and two-component SN la models have different assump- 
tions, they produce ICS Fe fractions in rough agreement. 
We do not explore more sophisticated theoretical mod- 
els that vary the delay-time distributions (DTD) of SN 
la (see Greggio 2005) to explain the Fe content in clus- 
ters (e.g. Maoz & Gal- Yam 2004; Ettori 2005) because 
of large uncertainties in the explosion mechanism and 
properties of SN la progenitors. 

Both the baseline and two-component SN la models 
predict that the ICS Fe contribution falls short of pro- 
ducing the observed Fe mass within the ICM. In some 
combination, we expect that additional enrichment chan- 
nels (such as metal-injecting winds from galaxies) and 
environment- and/or time-dependent SN la rates and 
yields will play roles in the final resolution of the question 
of ICM enrichment. We discuss metal loss from galaxies 
in §5.2, and focus here on whether our data indepen- 
dently support the claim of environment-dependent SN 
la rates (Mannucci et al. 2008). (Our sample lies at one 
redshift, so assessing time-dependencies from it is not 
possible.) 

One potential piece of evidence supporting the effect of 
environment is the trend between the model Fe fraction 
and a shown in Figure 9 (top panel). A Spearman rank 
test demonstrates that the null hypothesis (no trend) can 
be excluded with > 90% confidence. If the ICS contribute 
a fixed traction of the ICM Fe for all of our groups and 
clusters, then stars in more massive environments must 
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eject more Fc per stellar mass. If we attribute these 
Fc production differences entirely to differences in SN la 
rates, then we infer a factor of roughly three increase in 
that rate from the low-cr to high-cr limits of our range. 
This matches the factor of three increase between field 
and cluster SN la rates found by Mannucci et al. (2008). 
In this scenario, our baseline model overpredicts the Fe 
fraction coming from the ICS in group environments be- 
cause wc matched the present-day PEGASE SN la rate 
to that observed in clusters. In any case, the ICS alone 
cannot fully produce the ICM Fe. 

Given the variation among measurements of field SNe 
rates discussed above, and the likely dependencies on en- 
vironment and look-back time, the rates are a significant 
source of uncertainty in resolving the chemical enrich- 
ment problem in clusters. 

5.2. Metals from Galaxies 

The previous section indicates that the ICS is an im- 
portant, but not exclusive, contributor of Fe to the ICM. 
Therefore, we now estimate the contribution of galactic 
stars to the Fe ICM budget within r5oo for the seven 
Gonzalez et al. (2007) clusters for which we have optical 
photometry for all galaxies within rsoo. We first subtract 
the BCG+ICS contribution calculated using the baseline 
model from the measured ICM Fe value. The remainder 
of the metals must come from galactic stars. Once again 
we apply the baseline model (because the galaxy stel- 
lar population is dominated by old stars in early-types 
within rsoo) to predict how much Fe galaxies produce. 
We then calculate the galactic metal loss efficiency such 
that the sum of the weighted fractional contribution of 
the BCG+ICS and galaxies is 1. On average, the galax- 
ies need to lose 84^j^4% of the Fe that they produce for 
the sum of BCG+ICS and galaxies to account for 100% 
of the observed Fe (Figure 9; bottom panel). The choice 
of a single metal loss efficiency for all systems may not 
be appropriate as we over-predict the metals produced 
by two clusters, Abell 2984 and S 84. However, it is 
impossible to disentangle variations in SNIa rates (see 
above) and metal-loss efficiency with environment with 
the available data. 

This result does not consider that some fraction of all 
metals arc locked up in the current stellar populations of 
these galaxies. If the galaxies within r5oo have an average 
metallicity that is close to solar (see Portinari et al. 2004), 
approximately 25% of the metals must be locked up in 
stars. Therefore, given the uncertainties cited above, the 
baseline model satisfies two key constraints: 1) the com- 
bination of metals produced by the ICS and other cluster 
galaxies can account for the observed ICM Fe, and 2) the 
inferred metal loss fraction from galaxies (84ii[4%) does 
not confiict with the locked-up metal fraction (~25%). 

This baseline model requires that cluster galaxies lose 
a substantial fraction of the metals that they produce. 
While studies using sophisticated chemical models have 
shown that field early-type galaxies might lose ^ 80% 
of the metals (specifically Fe) they produce during their 
lifetime (Calura et al. 2006), other observational studies 
suggest a much lower metal loss fraction (measured from 
a elements) (~ 35%; Bouche et al. 2007, who do not 
distinguish between galaxies in different environments). 
This large discrepancy can at least be partly explained by 



Bouche et al. (2007) 's focus on a element measurements 
and their use of the instantaneous recycling approxima- 
tion (IRA) to compute galaxy metal loss eflSciency. While 
the IRA may off'er a reasonable approach for computing 
the production of a elements, it is not appropriate for 
computing the production of Fe because of the compar- 
atively longer timescale for Fe production (Rcnzini et al. 
1993). Fe is less likely to be locked-up in stars and more 
likely to be expelled into the IGM because the majority of 
Fe is produced after the epoch of star formation. There- 
fore, 35% metal loss fraction serves as a lower bound for 
the amount of Fe that can be lost. 

To examine the implications of the lower boimd of the 
metal loss fraction on our model, wc reverse the question 
and ask what the SN la rate needs to be if the Fe loss frac- 
tion from galaxies is only 35%. Given that condition, in- 
creasing the SN la rate to roughly its 2a upper- limit, i.e., 
by a factor of 1.8 x, succeeds in generating all the ICM 
Fe. Therefore, considering the large uncertainties in the 
SN la rates, we conclude that BCG+ICS+galaxy mod- 
els with galactic metal mass loss fractions ranging from 
^ 75% (the upper limit due to the locked-up fraction) 
down to 35% can reproduce the ICM Fe abundances. 

6. CONCLUSIONS 

We present X-ray surface brightness, temperature, and 
Fe abundance profiles out to ~ r5oo (~ 0.6 times r2oo, 
the virial radius) for a set of twelve nearby (z < 0.1) 
clusters with extensive optical photometry to quantify 
the relative contributions of intracluster stars (ICS) and 
cluster galaxies to intracluster medium (ICM) metal en- 
richment. We study systems with velocity dispersions of 
500 < (T < 1000 km s^^, measuring peak X-ray temper- 
atures and 0.5 — 2.0 keV luminosities within rsoo of 2 — 6 
keV and 0.3-3.4x10""^ ergs s-\ respectively. The T-a, 
Lx — <J, and Lx — T scaling relations for our systems are 
generally consistent with expectations. Two-dimensional 
X-ray surface brightness profiles reveal substructure in 
most of our clusters ranging from spiral-like structures 
to subclumps to central asymmetries. 

The majority (9 of 12) of the clusters, even those with 
isothermal profiles, show steep abundance gradients that 
approach the canonical abundance of ~ 0.3 Zpe.o at 
large radii. Four of these 9 have cool cores. The ra- 
dial variation in abundance makes necessary spatially- 
resolved spectroscopy to accurately quantify the radial 
distribution of metals, which is in turn necessary for self- 
consistent chemical evolution modeling. 

The stellar component consisting of the brightest clus- 
ter galaxy plus intracluster stars (BCG+ICS) contributes 
on average a non- negligible fraction, 31^9^%, of the 
ICM's Fe within rsoo for our baseline chemical evolution 
model. We also calculate the Fe yield for an alternate, 
two-component (prompt and delayed), empirical SN la 
model (Scannapieco & Bildsten 2005). Even though the 
baseline and two-component SN la models have differ- 
ent assumptions, they produce BCG+ICS average frac- 
tional Fe contributions that are similar: 3lt"% versus 

22+9% within rsoo, respectively. For the seven clusters 
in which we know the relative contribution of the BCG 
and ICS, the ICS contribute 80% on average of the com- 
bined BCG+ICS Fe, indicating that the ICS significantly 
enrich the intracluster medium within rsoo and must be 
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included in any enrichment model for the ICM. 

Because the BCG+ICS component cannot account for 
all of the Fe in the ICM, we consider the combined effect 
of BCG+ICS and other cluster galaxies on ICM enrich- 
ment. We then find that we can account for all the Fe 
within r5oo and that the required galactic metal loss frac- 
tion (841^4%) does not conflict with the fraction of met- 
als still locked up in galactic stars (~25%). While this 
metal loss efficiency is large, it is consistent with other 
estimates (^80%; Calura et al. 2006). It is also worth 
noting that, given the large uncertainties in the SN la 
rates, the required metal loss fraction might be signifi- 
cantly lower than this initial estimate. For example, if we 
increase the rate of SN la in cluster galaxies by a factor 
of 1.8, to its upper 2<t bound (Mannucci et al. 2008), we 
can produce all the ICM's Fe with a galactic metal loss 
fraction of only ^35%, which is consistent with other 
metal loss fraction estimates (Bouche et al. 2007). As 
a result, it is possible to make a full accounting of the 
ICM's Fe — without resorting to extreme assumptions 
about the stellar initial mass function or pre-enrichment 
— for a large range of plausible galactic metal mass loss 
fractions (35 - 75%, the upper limit due to the locked- 
up fraction). Reducing the uncertainties in the SN la 
rates, including measurements of the likely dependencies 
on redshift and galaxy environment, is the critical next 
step in understanding the enrichment history of the in- 
tracluster medium. 
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TABLE 1 

Observed X-ray clusters and their properties 



Name 


(^) 


BM 


O" 

(km s-l) 


rsoo 
(Mpc) 


r200 
(Mpc) 


(1044 


Lx 

ergs s 1) 


Tpeak 

(keV) 


Abell 
Abell 


496 
1651* 


0.0329 
0.0845 


I 

I-II 


743+4?, 
990+"" 


1.033 
1.397 


1, 
2 


.599 
.164 


1 

3, 


01+0.06 
■"-'--0.08 
Q7+0.16 
—0.16 


4.39 
6.16 


± 
± 


0.08 
0.22 


Abell 


2811* 


0.1079 


I-II 


86ot;s;; 


1.202 


1 


.862 


2, 




5.76 


± 


0.20 


Abell 


2877 


0.0247 


I 


999+71 


1.406 


2, 


.177 






3.27 


± 


0.14 


Abell 


2984* 


0.1042 


I 


490+112 


0.658 


1, 


.020 





QQ+0.02 

■'^"^-0.03 


2.01 


± 


0.05 


Abell 


3112* 


0.0750 


I 


940tl« 


1.324 


2 


.050 


2 


7Q+0.04 
• '^-0.04 


5.24 


± 


0.11 


Abell 
Abell 


3693* 
3705* 


0.1237 
0.0906 


III 


losoiiis; 


1.452 
1.426 


2 
2 


.249 
.209 





on+0.04 
•"■^-O.OT 


4.21 
3.71 


± 
± 


0.14 
0.22 


Abell 


4010* 


0.0963 


I-II 




0.867 


1, 


.343 


1, 


46+0-08 
■^"-0.14 


4.44 


± 


0.34 


Abell 


4059 


0.0475 


I 


6531^4 


0.901 


1, 


.396 


1, 


Qo+o.oa 

■"^"^-0.07 


4.14 


± 


0.13 


Abell 


S 84* 


0.1100 


I 


520tl«° 


0.705 


1 


.091 





7f;+0.13 
• '°-0.08 


4.28 


± 


0.18 


Abell 


S540 


0.0358 


I 


760±* 


1.057 


1 


.638 


0, 


qo+0.01 

•"^•^-o.os 


2.65 


± 


0.08 



Note. — BM = Bautz-Morgan type. Velocity dispersions are from Zaritsky et al. 
(2006) or calculated for this paper using the method outlined in Zaritsky ct al. (2006) using 
radial velocities from the NASA Extragalactic Database. Luminosities arc measured in the 
0.5 — 2.0 keV range out to r2oo. Our data do not extend that far in radius and therefore we 
extrapolate our surface brightness fits to obtain the enclosed X-ray luminosity. The quoted 
temperatures are the peak values obtained from our spectral fits. We choose the peak value 
as it best reflects the depth of the cluster potential, especially for clusters with a cool-core. 

* Gonzalez et al. (2005) clusters. 



TABLE 2 

X-RAY DATA PROPERTIES 



Name 


Obs IDi 


Date 


Type2 


Filter 


tMOS 

(ks) 


tMOS.fl 
(ks) 


router 

(arcmin) 


Abell 


496 


0135120201 


2001 Feb 01 


XSA 


Thin 


29.5 


16.1 


10 


Abell 


1651 


0203020101 


2004 Jul 01 


GO 


Thin 


15.2 


8.0 


6 


Abell 


2811 


0404520101 


2006 Nov 28 


GO 


Medium 


24.4 


22.7 


5 


Abell 


2877 


0204540201 


2004 Nov 23 


XSA 


Thin 


21.9 


20.1 


10 


Abell 


2984 


0201900601 


2004 Dec 27 


XSA 


Thin 


29.2 


27.2 


4 


Abell 


3112 


0105660101 


2000 Dec 24 


XSA 


Medium 


23.3 


22.4 


6 


Abell 


3693 


0404520201 


2006 Oct 14 


GO 


Medium 


34.6 


29.4 


3 


Abell 


3705 


0203020201 


2004 Oct 06 


GO 


Medium 


24.4 


13.9 


6 


Abell 


4010 


0404520501 


2006 Nov 13 


GO 


Thin 


19.6 


18.5 


6 


Abell 


4059 


0109950101 


2000 Nov 24 


XSA 


Thin 


29.4 


12.4 


9 


Abell 


S 84 


0201900401 


2004 Dec 04 


XSA 


Thin 


34.6 


17.7 


4 


Abell 


S540 


0149420101 


2002 Oct 11 


XSA 


Medium 


18.0 


11.1 


8 



iXMM-Newton observation identification. 

^XSA denotes data from the XMM-Newton Science Archive; GO denotes data from our 
Guest Observer programs. 



The Enrichment of the ICM 



TABLE 3 
Best fit 2D beta model parameters 



Name 




rx 












(kpc) 


(cm-3) 




Abell 


0496 


473+S SS^ 


21 1+°? 


033+2 SS? 

U.UOO_Q QQ2 


6.6 


Abell 


1651 


549+S SS2 


90 7+2 5 
' —2.5 


oio+2S2^ 


1.3 


Abell 


2811 




141.713-7 


006+° °°3 

U.UUD_g QQ3 


1.1 


Abell 


2877 


0.28 


0.38 




2.9 


Abell 


2984 


U.OUi_,J QQg 


33.612-2 


012+° °°l 

U.UlZ_g ggj^ 


1.2 


Abell 


3112 


564+"'""-'* 

U.dO'i_„ (,(,3 


51.6li;i 


020+° °°l 

U.UZU_Q ,„)j 


1.7 


Abell 


3693 




82-011;; 


00^2+°-'"""* 


1.7 


Abell 


3705 










Abell 


4010 


U.OUU_„ (,(,3 


24 2+"-'^ 


031+0-""2 

U.U01_g gg2 


2.8 


Abell 


4059 


497+0. 004 


43 1+11 




1.1 


Abell 


S 84 




66.8l|;5 


0078+° °°°'^ 


1.3 


Abell 


S540 


49+0 01 


39.4l^-f 


0074+0-'""'8 

U.UU(*_g QQQ7 


1.2 



1 One-dimensional value normalized by the degrees of freedom 
(DOF). 



TABLE 4 

X-RAY DERIVED CLUSTER HYDROGEN GAS MASS 



Name 


Tspec 


rspec/rBOO '^-H.i-spec 


(1013 Mq) 




(kpc) 


(1012 Mq) 



Abell 


0496 


392. 


,7 


0, 


.38 


6.112:1 


2, 


,851°; 


Abell 


1651 


570. 


.5 


0. 


.41 




6, 


,601°; 


Abell 


2811 


590. 


.5 


0. 


.49 


15 5+°-2 

ld.d_g 2 


3, 


,901°; 


Abell 


2984 


458. 


,9 


0, 


.70 


4.4lg| 


0, 


,771°; 


Abell 


3112 


511, 


,9 


0, 


,39 




4, 


,451°; 


Abell 


3693 


399. 


,7 


0, 


,28 


5.ol?,:l 


3, 


■651°; 


Abell 


4010 


642. 


.0 


0. 


,74 


12.21;- 


1, 


,931"; 


Abell 


4059 


501. 


,9 


0, 


,56 


8.4l°;l 


2, 


,091°; 


Abell 


S 84 


481. 


,0 


0, 


,68 


6 0+°-^ 

D.U_g 2 


1, 


,011°; 


Abell 


S540 


339. 


,8 


0, 


,32 


n QQ+O.OS 

Z.OO_g gg 


1, 


•4ll°o: 



1.14 



TABLE 5 

Cluster gas mass comparison with literature 



Name 


rsoo 
(kpc) 


Ref M9,^5„„ 

(1013 Mq) 


Our Mg.rsoo 
(1013 Mo) 


Ref 


Cosmology 


Abell 496 


1420±30 


6.9 


8-4lg;! 


1 


SCDM50 




12401^° 


6.8l°g;^ 


6.8l°g;l 


2 


SCDM50 




668±96 


1.2±0.3 


i-9l°o;! 


3 


LCDM70 


Abell 1651 


1550±50 


12.4 


14.01°:! 


1 


SCDM50 




17301^° 


15-illi 


16.4+0-9 


2 


SCDM50 


Abell 3112 


1440l?°g 


9.4 


9 4+it 


1 


SCDM50 




1530111° 


I0.5ll;3 


10.2lJ;| 


2 


SCDM50 




1048±317 


2.8±0.6 


4 3+21 


3 


LCDM70 


Abell 4059 


1400l«g° 


6 4+°-*' 


7-4l°g;i 


2 


SCDM50 




1247±304 


2.4±0.6 


4.611:? 


3 


LDCM70 


Abell S540 


108011^° 


2 4+°* 
^■^-0.6 


2.6l°g:| 


2 


SCDM50 



References. — (1) Castillo-Morales & Schindler (2003); (2) Reiprich 
(2001); (3) PifTaretti et al. (2005) 

Note. — rsoo values taJjen from references. 
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TABLE 6 

DePROJBCTED BCG and ICS Fe contribution out to rap, 



Name MFe,r,p„c '!^Fe,BCG,ispec ^Fe,IC S .r^pec ^Fe,BCG+ IC S .i^p^c '^Fe,BCG+IC S /^^^^ Fe 

+0.4 _ _ ^ o 1 n2+°-i^ 

-0.4 J^-'J^-O.IT 

7.5 3.5 11.0 0.651°;^^ 

ytl;l 1-8 8.5 10.2 0.68lo;ig 

1.1 8.9 10.0 2.03t°;4g 

ll'l 1.0 10.6 11.6 0.96l„;^7 

t\:l 1-7 4.9 6.7 l-48lo:43 

+il 1.7 8.3 9.9 0.60l[J-J^ 

6.4 0.861°;^^ 

7.5lJ-^ 0.4 7.7 8.1 l-08lo25 

2 4+°-^ - - 6 6 2 77+0'S4 

^■^-0.3 ^•"-0.54 



Abell 


0496 


Abell 


1651 


Abell 


2811 


Abell 


2984 


Abell 


3112 


Abell 


3693 


Abell 


4010 


Abell 


4059 


Abell 


S 84 


Abell 


S540 



4.9 
12.1 
4.5 
16.6 
7.4 



Note. — All Fc masses are in 10^ Mq units. These results are derived from our baseline PEGASE model (Fioc 
& Rocea-Volmerange 1997), which reproduces the present-day current cluster early- type SNIa rate (Mannucci 
et al. 2008). For the two-component SN la model values, divide the BCG and ICS Fe masses, and the Fe 
fractional contributions by 1.39. 



TABLE 7 

Deprojected BCG and ICS Fe contribution out to rsoo 



Name 


^Fe,i5oo 


MFe,SCG,r5oo 


MFe,/CS,r5oo 


MFe,SCG+/CS,r5oo 


^Fe,BCG+IC 


Abell 0496 


no + 3.4 
ZO.O_3 g 






6.5 


0.29: 


Abell 1651 




7.7 


4.5 


12.2 


o.2o: 


Abell 2811 


39.0t|:l 


1.8 


9.5 


11.3 


0.29: 


Abell 2984 


5 q+lO 
'^■-'-1.0 


1.1 


10.0 


11.1 


1.88: 


Abell 3112 


36.0+g 7 


1.0 


12.1 


13.1 


0.36: 


Abell 3693 


34.6tl^:«3 


1.8 


8.1 


9.8 


0.28: 


Abell 4010 


24.9lt^ 


1.7 


8.9 


10.5 


0.42: 


Abell 4059 


14.7ti:| 






6.9 


0.47: 


Abell S 84 


12.3i|2 


0.4 


8.1 


8.4 


0.69: 


Abell S540 


c 7+2.3 
^•'-1.5 






7.3 


0.55: 



+0.06 
-0.06 
+0.06 
-0.05 
+0.08 
-0.08 
+0.43 
-0.43 
+0.09 
-0.09 
+0.12 
-0.11 
+0.11 
-0.11 
+0.14 
-0.14 
+0.17 
-0.16 
+0.15 
-0.15 



Note. — All Fe masses are in 10® Mq units. These results are derived from our baseline PEGASE model 
(Fioc & Rocca-Volmcrange 1997), which reproduces the present-day current cluster early- type SNIa rate 
(Mannucci et al. 2008). For the two-component SN la model values, divide the BCG and ICS Fe masses, 
and the Fe fractional contributions by 1.39. 
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A2S1 1 



I 



-lil 



IJ.J 




BBSS BB.ta 
(dag) 




1Sa.DD 194.^ 
RA. (dag) 




ia.E5 iD.aa 



-418 

'W -15.9 

1 




A31 12 



17,7a 17.55 17.40 



17,a 



40.2 f "j:-... 



I-' - .."Ih: 



i.-V"r. " a: 



■I'-.v.- 



33j[)D 3Z.B& S.TQ 




4a.es 49.50 43.35 



-34^ 



1 



-34J : 



308.70 3Da.5S 




11* C^) 



510vS 



A4<910 



3 -3Bi ^. 




5aZ95 SaZfiD 
(d^ 



-34,e 



-34,7 



* -34.sif.;^ 




351,29 aD,ia 



AS 84 



AS540 



».5 



-29.7 




1£4a 12,30 



KOa B4,fl0 



Fk;. 1. 15' by 15' Digital Sky Survey (DSS) images with overlaid 0.3 — 8.0 keV X-ray contours (particle background subtracted) for 
all observed clusters. The DSS images are stretched by first carrying out a 3cr from the mean clip of all pixels. The lower and upper limits 
of the stretched image are chosen to be the mean of the clipped pixels and the mean plus 3a, respectively. The point sources have been 
removed in the smoothed X-ray maps used to generate contours. Contour intervals are logarithmically spaced, spanning 5 intervals from 
the lowest value of 2x10^'' counts s^^ arcscc^'^ to the peak value in the smoothed map. In general, the clusters arc "relaxed," i.e. arc 
roughly spherical, and the contours arc centered about the brightest cluster galaxy (BCG). Abell 2877 and 3705 are exceptions. The X-ray 
centroid of the diffuse emission of Abell 2877 is not centered on the BCG. Abell 3705, which is the only BM Type III cluster in our sample, 
is not spherical and may be in the process of a merger with the NW clump. 
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200 400 600 800 1000 200 400 600 800 200 400 600 800 
Projected Rodius (kpc) Projected Rodlus (kpc) Projected Rodius (kpc) 



Fig. 2. — Surface brightness [0.3 — 8.0 keV] profiles and ID fit residual in units of x (left), temperature profiles (middle), and metallicity 

profiles (right) of the sample of clusters. The particle background has been subtracted from the surface brightness profiles, leaving only 
the cosmic background. Several clusters with cool cores have temperature profiles that rise from the center, peak, and then drop at larger 
radii. A majority of the clusters have radial abundance gradients that flatten to Z ~ 0.3 Zp^ q at large radii. The two clusters (Abell 1651 
and 2811) with no obvious abundance gradients are also isothermal. Abell 3705 has such low signal-to-noise that its surface brightness flt 
fails to converge, as shown by its residual plot. 
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Fig. 2. — Continued. 




Fig. 3. — Abell 3112 2D surface brightness profile (10' on each side) fitting residuals. The data arc smoothed by a 25" gaussian kernel 
to maJce the residuals easier to see. The color bar is in units of counts. The original data with the point sources removed are shown in the 
top left. This plot illustrates the ability of different types of models to fit the observed surface brightness. Spherically symmetric model 
fits (top right) result in asymmetric residuals because of the elliptical nature of cluster emission; a single component elliptical beta model 
(bottom left) fits fairly well, but only marginally better than the spherical fit. Adding an additional beta model (bottom right) to the fit 
again only improves the fit marginally. 
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Fig. 4. — 0.3 — 8.0 kcV 2D surface brightness fit residuals (15' on each side) for observed clusters. The data are siiiootlied by a 25" 
gaussian kernel to make the residuals easier to see. The color bar is in units of counts. These residuals are not normalized by the total 
observation time. Several cluster residuals have a distinct quadrupole pattern that is an artifact of carrying out a circularly symmetric 
model fit to an intrinsically elliptical surface brightness distribution. However, several clusters reveal real minor substructure that is clearly 
asymmetric. An extreme case is Abell 496, where there is an apparent spiral structure leading to the center of the cluster. AbcU 3112 
shows an extended tail to the southwest. AbcU 3705's residuals show a very poor fit due to the cluster's high asymmetry and the data's low 
signal-to-noise. In Abell 3693, there is a bright clump to the southeast. AbcU 4059 has a filamentary structure near the center. Overall, 
the observed substructure does not significantly impact our analyses, as evidenced by the acceptable values we obtain for most cluster 
surface brightness fits. 
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Fig. 5. — Testing systematic errors arising from the spherieal symmetry assumption. In the left panel we show the azimuthally averaged 
surface brightness profile for a range of shapes (the legend provides tx,ty,tz for each model) . Although the different clusters are offset 
in surface brightness the profiles are similar in shape. The right panel demonstrates this to be the case after we normalize the profiles at 
r = 25. These two plots suggest that intrinsic shape differences should not lead to dramatic differences in profile parameters other than the 
normalization, except at very large or small radii. The bias is for elliptical clusters to be fainter than the corresponding spherical cluster. 
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Fig. 6. — Left: Differences in the 0.1 — 2.4 keV luminosities determined by this work (Lx) and those from REFLEX {Lueflex', 
Bohringer et al. 2004) within an aperture defined by the REFLEX ROSAT observations. The measurement differences for most clusters 
fall within the 20% envelope, although our measurements tend to be brighter than theirs at higher velocity dispersions. Right: Differences 
in temperatures determined by this work (Tx) and those from the BAX database (Tg^x) for seven clusters. Three out of four clusters 
(Abell 496, 3112, and 4059) with systematically lower temperatures in BAX have cool-cores. We expect the BAX temperatures, typically 
derived using emission-weighted methods, to be systematically lower than our peak temperatures in such cases. However, according to our 
spectral fits, the fourth cluster, Abell 2811, is isothermal, suggesting its discrepancy is real. 
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400 600 800 1000 1200 

a (km s~') 

Fig. 7. — Gas mass within rsoo plotted against cluster velocity dispersion. The diamond points are our clusters and the triangular points 
are the seven Vikhlinin et al. (2006) X-ray clusters with known velocity dispersions. The dashed line is a power law fit to the Vikhlinin et 
al. (2006) data, while the solid black line is the fit to our data. The fits agree closely. We caution when interpreting results from the lowest 
velocity dispersion clusters, because the relationship used to estimate the rsoo values is poorly constrained in that velocity dispersion range. 
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Fig. 8. — Left: Cluster temperature (derived from peak cluster temperature) versus cluster velocity dispersion. The text beside each 
datapoint is the cluster's name. The solid line is the T-cr fit taken from Wu et al. (1999); the dotted lines represent Itr errors of the fit. The 
majority of the clusters lie within the \a errors of the fit. The arrows point to alternate velocity dispersions for Abell 2877 (after removing 
its lower velocity dispersion peak) and Abell S540 (after clipping its velocity wings). Middle: 0.5 — 2.0 kcV X-ray luminosity plotted against 
velocity dispersion. The solid line is the best fit excluding the two outliers, Abell 3693 and S540. Right: 0.5 — 2.0 kcV X-ray luminosity 
plotted against velocity dispersion. The solid line is the best fit for all clusters. These X-ray-optical relations further establish the validity 
of the X-ray reduction and analysis techniques employed here. 
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Fig. 9. — Top: Fraction of the ICM's Fe within rsoo contributed by the BCG+ICS component as a function of velocity dispersion. The 
dotted line at 31% is the weighted average of the fractional Fe contribution of the full sample. Averages are computed using l/error^ as 
weights. The Fe fractional contribution changes significantly with velocity dispersion (see text). However, the rsoo values for two lowest 
velocity dispersion clusters (Abell 2984 and Abell S 84) are only constrained by a single low velocity dispersion point in the calibration of 
the rsoo — c relation (Gonzalez et al. 2007), and wc suggest caution in interpreting results from these two data points. Middle: Fraction 
of the ICM's Fe within rsoo contributed by non-BCG cluster galaxies as a function of velocity dispersion. We assume 84% metal loss, the 
value required to generate on average all the ICM's Fe (sec bottom panel) from the combined eoiitiibutiou of the BCG+ICS and galactic 
stellar components. We plot only the subset of clusters for which the galactic stellar component is well measured (Gonzalez et al. 2007). 
The dotted line at 69% is the weighted average of the galactic contribution. The galaxies alone cannot account for all the metals in the 
ICM even if they have a 100% metal mass loss efhcieney. Bottom: Fraction of the ICM's Fe within rsoo contributed by all the stars — in 
and out of galaxies — as a function of velocity dispersion. Only clusters with both BCG+ICS and galactic light measurements are plotted. 
A galaxy metal loss fraction of 84% is chosen so that the weighted average of the total stellar mass contribution (dotted line) is 1 , thereby 
accounting for all the measured Fe. While this metal loss fraction is large, it is consistent with some estimates (~80%; Calura et al. 2006) 
and with the fraction of metals not still locked up in galactic stars (~75%). If we increase the rate of SN la in cluster galaxies by a factor 
of 1.8, to its upper 2a bound (Mannucci et al. 2008), we can produce all the ICM's Fe with a galactic metal loss fraction of only ^35%, 
which is consistent with other metal loss fraction estimates (Bouche et al. 2007). 
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APPENDIX 
NOTES ON INDIVIDUAL CLUSTERS 
Abell 496 

This is a nearby high signal-to-noisc, mildly elliptical, cool-core cluster with emission that overfills the XMM-Newton 
field-of-view. This complicates surface brightness analysis because the flat cosmic background is masked by cluster 
emission that is at least an order of magnitude stronger within the field-of-view. The initial surface brightness fits 
converge at zero background due to the combination of asymmetric structure near the core and bright cluster 
emission compared to the background. In Figure 4 we show the spiral-like structure that remains in the central parts 
of the cluster after subtracting off the best fit beta model. This feature is also present in Chandra exposures of the 
cluster core. It is suggestive of cooler gas spiraling into the cluster's center, but this can also be explained by chance 

alignment of higher density gas clumps. We correct the poorly fit background by carrying out a new fit with the 
central 4' excised. This removes the central asymmetries, resulting in a fit that converges on a non-zero background. 

We then fix the background in the original surface brightness fit to the newly derived background values and fit 
again for the beta model parameters. We use these values in our further analyses. We also excise the central 20" of 
this cluster to remove the brightest cluster galaxy and fit only the extended cluster component of the beta model. 
Despite our attempts at improving the fit, the surface brightness residual, see Figure 2, is large, showing a systematic 
trend with radius. A possible explanation for this is the fainter asymmetric structure observed in the residual north 

of the cluster center. 

Abell 1651 

One of the higher redshift and most massive clusters in our sample, this is the best-behaved cluster with reasonable 
signal-to-noise, even though half of the exposure time was ruined by flaring events. The X-ray contours are circular 
and are centered on the brightest cluster galaxy. The cluster temperature proflle is isothermal and, as expected, its 
metallicity proflle shows no gradient. The surface brightness fit is consistent with the data and shows no systematic 
variations with radius. Two-dimensional fit residual in Figure 4 shows a mild asymmetric structure centered around 

the cluster core, which is an artifact of the circularly symmetric models as the cluster is slightly elliptical. The 
relative contributions of the BCC and ICS to the BCG-I-ICS light is poorly understood for this cluster. This is the 
only cluster in our sample that can be well fit by a single de Vaucouleurs profile (Gonzalez et al. 2000), making the 
distinctions between BCG and ICS somewhat arbitrary in the two component de Vaucouleurs fit (Gonzalez et al. 

2005). 

Ahell 2811 

This is another well-behaved cluster whose surface brightness profile is close to circiilarly symmetric. However, the 
difl[use cluster emission is offset from the brightest cluster galaxy by 27" in the northwest direction. There are no 
systematic trends observed in the surface brightness residuals. Just like Abell 1651, this cluster's temperature profile 
is isothermal and exhibits no metallicity gradients. This cluster is a significant outlier in the Lx — T relation and 
appears to be underluminous for its temperature. When we compare to the BAX results (Sadat et al. 2004), Figure 
6, we see that our derived cluster temperature is higher than the BAX temperature. We also notice that this cluster 
is hotter than what is expected for its luminosity in our Lx — T relation, suggesting that these discrepancies are real. 

Ahell 2811 

The lowest redshift cluster in this sample, Abell 2877 is one of the least well-behaved clusters. Its surface brightness 
profile has an uncharacteristically low (3 value of 0.27, and it is not well-fit by a beta profile. The X-ray image shows 
a very luminous central source surrounded by much fainter extended diffuse emission. The central source, which is 
coincident with the brightest cluster galaxy, is offset south from the center of the diffuse emission. We excise a 47" 
radius region around the central source to remove all of its light and fit just the extended emission. Despite our 
attempts to remove any sources of problems for our fit, the fit consistently converges to unphysically small core 
radius and small beta value. A NASA Extragalactic Database^ (NED) query for this cluster shows that there are two 
superimposed clusters separated by lS.z ~ 0.006 in redshift space, which is a likely explanation for the poor surface 
brightness profile fits. A velocity histogram of the member galaxies unsurprisingly shows a second, albeit smaller, 
peak, likely leading to an overprediction of the cluster velocity dispersion. This cluster is an outlier in the T — a 
relation, suggesting that cr is indeed overestimated. Removing the smaller velocity peak from the distribution and 
recalculating a moves A2877 to within the errors of the relation. We also center the spectral extraction radii used for 

our spectroscopic fits about the the extended emission. This cluster exhibits a metallicity gradient, and the 
abundance falls precipitously at large radii, which is not observed in the other clusters. For these reasons, we remove 

this anomalous, apparently unrelaxed cluster from our analysis. 

Abell 2984 

This is the lowest velocity dispersion cluster in the sample and also the smallest in physical extent (see Figure 1). The 
cluster is not axis-symmetric, though the surface brightness fit does a good overall job fitting the cluster, as indicated 

® http://nedwww.ipac.caltech.edu/ 
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by the fit residuals. This cluster does not look peculiar in the various T — Lx ~ cr relations. Its gas mass also appears 
consistent with that expected from the relation derived from gas mass data from Vikhlinin ct al. (2006), Figure 7. 
However, the BCG+ICS makes an unphysical fractional Fe contribution (greater than unity) to the ICM metal 
budget (see §5.1). The origin of this overprediction is unclear. Figure 6 suggests that we underpredict the luminosity 

of this cluster, which will in turn lead to an underprediction of gas mass and metal mass. It may also be that the 
intracluster light measurements are flawed and conspire to produce too large a contribution to the ICM metals, or it 
could be that our scaling of the PEGASE model to match the cluster early-type SN la rate may not be applicable for 
this large group. If we neglect this rescaling, then the BCG+ICS contribution falls to a physically believable value 
(less than unity). Note that this cluster has both a metallicity gradient and an isothermal temperature profile. 

Ahell 3112 

This cluster is relatively well-behaved, though it is one of the more elliptical clusters in the sample. It has some 
extended emission from what appears an infalling clump south of the cluster center. This clump, however, has 
negligible impact on the overall cluster surface brightness distribution and is not evident in the X-ray contour plots. 

Overall, the beta profile does a good job of describing the cluster surface brightness profile. There is some 
discrepancy between the fit and the data at small radii, but that is expected when fitting a circularly symmetric 
model to an elliptical distribution. We do excise the inner 35" to remove the cool-core region from the fit. This 

cluster does not appear anomalous in all of our diagnostics. 

Ahell 3693 

This is the highest redshift cluster in our sample and has moderate signal-to-noisc data. Even though the diffuse 
cluster emission appears circular, we observe a large clump to the southeast of cluster. This structure is evident in 
the surface brightness profile and fit residuals at a radius of 900 kpc from the cluster center. The velocity histogram 

has multiple peaks, suggesting that the velocity dispersion of this cluster may be overestimated. This cluster is a 
significant outlier in both the T — a and Lx — c relations, which is no surprise, but it is also somewhat discrepant in 
the Lx — T plot. One possible explanation is that this system is undergoing a minor merger, which is not likely to 

significantly impact our analyses. 

Ahell 3105 

This is the least spherically symmetric and lowest signal-to-noise cluster in our sample, making it impossible for a 
surface brightness fit to converge. In the X-ray contour plot, two chimps of equivalent size dominate the X-ray 
emission. For a simple circularly symmetric beta profile, the fit parameters converge to senseless values and produce 
very large fit residuals. Multiple approaches, including excising a large region around the smaller clump to the 
northwest of the cluster center, fail. For this reason, we remove the cluster from further analysis. However, it is 
instructive to look at the various diagnostic plots, because we are able to extract a few X-ray parameters from the 
data. This cluster is an outlier in the T — a plot. Its velocity histogram has multiple peaks hinting at substructure 
and suggesting that our a is an overestimate. This unrelaxed system is the sole Bautz-Morgan Type HI cluster in 

our sample. 

Ahell 4010 

The X-ray data for this cluster is of moderate signal-to- noise. The surface brightness map looks regular with the 
exception of a slight asymmetry to the southeast in the last X-ray contour in Figure 1 and in the surface brightness 
fit residual 300 kpc from the cluster center. We excise the inner 35" to remove the cool-core and fit the extended 
X-ray component. The beta model fit does a reasonable job in fitting the surface brightness profile. The cluster 
appears "normal" in the T — a, Lx — cr, and Lx — T diagnostics. 

Ahell 4059 

We obtained high signal-to-noise X-ray observations of this cluster. The cluster surface brightness distribution 
appears smooth and is mildly elliptical. The beta model fit shows no systematic trends in the residuals, and the 
cluster falls within all of the trends found in our diagnostic plots. There is nothing to suggest this cluster is peculiar. 
We do not excise the central region of the cluster because our beta profile did not show any large residuals in the 

cluster center. 



Ahell S 84 

Our moderate signal-to-noise data show this cluster to be regular and circular. Its surface brightness fit is 
well-behaved. There is a hint of substructure 200 kpc southeast of the cluster center, which appears as a small bump 
in the residuals. This cluster is roughly consistent with the T — a relation, but is underluminous in the Lx — T 

plane. One possible explanation is that the temperature is overestimated for this cluster. Even though the 
temperature profile is isothermal, the metallicity profile has a steep gradient. 

Abell S540 

The relatively high signal-to-noise data show the highly elliptical nature of this poor cluster. We excise the central 
35" to remove the cool-core before conducting our surface brightness fit. There is no systematic trend in the fit 
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residuals. Despite the high velocity dispersion, the cluster temperature suggests that this system is a rich group or 
poor cluster. It is a 3(t outlier to the T — a and Lx — c fits. The cluster, however, falls on the mean Lx — T trend 
that we observe for our sample. The velocity histogram for the member galaxies is poorly-sampled, so it is likely that 

the velocity dispersion is overestimated. 



